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ABSTRACT 


Analysis  of  water  column  stability  was  performed  using  shipboard  and  submarine 
shear  and  density  profile  data  acquired  by  U.S.S.  DOLPHIN  and  R  V  ACANTA  in 
October.  19S4  in  the  vicinity  of  Monterey  Bay,  California.  Data  was  acquired  utilizing 
CTD  and  acoustic  doppler  profiler  (ADVP)  instruments.  The  upper  ocean 
thermohaline  structure  and  water  column  stability,  over  a  10km  square  domain,  was 
determined  from  repeated  measurements  of  the  conductivity,  temperature,  and  velocity 
to  a  depth  of  115m.  The  temporal  and  spatial  variation  in  the  analyzed  fields  of 
temperature,  salinity,  density,  and  velocity  are  compared  with  the  constructed  profiles 
of  the  static  stability  parameter  (E)  and  gradient  Richardson  number.  The  analyzed 
fields  were  in  turn  compared  with  the  larger  scale  forcing  factors  of  coastal  upwelling, 
current  systems,  bottom  topography,  and  internal  waves.  The  stability  the  water 
column  was  found  to  be  both  statically  and  dynamically  stable  with  the  exception  of 
thin  patches  of  instability  which  were  determined  to  be  the  result  of  double  diffusive 
processes. 
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I.  INTRODUCTION 


A.  OBJECTIVE 

The  objective  of  this  thesis  is  to  examine  the  static  and  dynamic  stability  of  the 
upper  115m  of  the  water  column  in  terms  of  the  spatial  and  temporal  variations  in  the 
density  and  velocity  fields.  Observations  are  based  upon  profile  data  acquired  by  the 
L'.S.  Navy  research  submarine  L'.S.S.  DOLPHIN  and  the  Navy  Postgraduate  School 
research  vessel  ACANIA  during  a  series  of  measurements  conducted  in  Monterey  Bay 
in  the  vicinity  of  the  Monterey  Canyon.  From  these  measurements  vertical  profiles  of 
the  static  stability  parameter  and  gradient  Richardson  number  were  constructed. 
Comparison  of  these  profiles  with  the  temperature,  salinity,  density  and  velocity  fields, 
as  well  as  the  synoptic  meteorological  analysis  and  local  bathymetry  were  made  in 
order  to  determine  the  effects  of  atmospheric  and  oceanographic  forcing  on  the 
stability  of  the  upper  100m  of  the  water  column. 

The  purpose  for  the  data  collection  was  an  experiment  lead  by  T.R.  Osborn 
(Osborn  and  Lueck,  1985)  for  submarine  measurement  of  turbulence  microstructure  in 
relation  to  vertical  gradients  of  density  and  velocity.  In  order  to  provide  a  larger  scale 
context  to  the  submarine  measurements,  concurrent  CTD  and  acoustic  doppler  velocity 
profiler  (ADVP)  measurements  were  conducted  by  the  R/V  ACANIA.  Measurements 
were  acquired  within  a  10km  'box’  centered  at  36°  44.5'N,  122°  03.0AV  over  a  35  hour 
period  from  the  afternoon  of  3  October  to  the  morning  of  5  October,  1984  (Fig.  1.1). 
During  this  period  R/V  ACANIA  completed  seven  circuits  of  the  edge  of  the  survey 
area  while  U.S.S.  DOLPHIN  conducted  two  dives  within  the  interior  of  the  survey- 
area.  The  separation  between  the  two  research  vessels  was  to  ensure  a  margin  of  safety 
for  the  submarine  in  the  event  of  emergency  surfacing.  In  order  to  maximize  the 
horizontal  resolution  of  the  temperature,  salinity  and  velocity  measurements,  a  CTD 
was  continuously  profiled  while  concurrently  recording  velocity  data  from  a  shipboard 
ADVP.  The  resulting  space,  time  series  allowed  the  vertical  structure  of  the  upper 
ocean  to  be  determined  over  scales  of  0.5m  to  100m  with  a  near  horizontal  sampling  of 
500m.  Data  collection  and  processing  are  discussed  further  in  Chapter  3. 
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B.  BACKGROUND 


The  waters  of  the  central  California  coast  are  subject  to  many  types  of  dynamic 
forcing  which  result  in  spatial  and  temporal  variations  in  the  thermohaline  structure 
and  velocity  fields  which  in  turn  affect  the  stability  of  the  water  column.  The  local 
forcing  is  dominated  by  the  eastern  boundary  current  system  (California  current  and 
countercurrent),  atmospherically  forced  coastal  upwelling,  tidal  action,  and  bottom 
topography. 

The  California  Current  system  is  composed  of  a  broad,  weak  equatorward  flow 
near  the  surface  (California  Current)  and  a  narrower,  submerged  poleward  counterflow 
(California  Countercurrent)  adjacent  to  the  continental  slope  (Wickham,  ei.  al.  19S6). 
In  the  period  from  October  to  March  there  is  typically  a  decline  in  the  intensity  of  the 
northwesterly  wind  flow  near  the  coast.  This  acts  to  reduce  coastal  upwelling  thereby 
allowing  the  countercurrent  to  surface  whereupon  it  is  known  as  the  Davidson  Current. 
During  this  time  period,  the  Davidson  Current  is  present  at  all  depths  shoreward  of  the 
California  Current.  It  is  characterized  by  a  relatively  thick  homogeneous  upper  layer 
and  poleward  flow  of  30  cm  s  or  less  (Blumberg,  1975). 

Coastal  upwelling  intensifies  in  the  spring  with  the  establishment  of  a  strong 
atmospheric  high  pressure  system/ridge  off  of  the  southern  and  central  California  coast 
and  the  concurrent  development  of  a  thermally  induced  low  pressure  system  (heat  low) 
over  the  interior  region  of  the  state.  The  northeast-southwest  orientation  of  the 
pressure  gradient  results  in  fairly  consistent  northwesterly  winds  along  the  coast 
through  the  spring  and  summer  months.  This  onshore  flow  causes  offshore  Ekman 
transport  in  the  ocean  resulting  in  coastal  upwelling.  The  upwelling  process  is  very- 
sensitive  to  the  atmospheric  forcing  and  thus  varies  directly  in  intensity  with  the 
surface  wind  field  with  variations  observed  on  the  order  of  several  days  or  less 
(Breaker.  1983).  Upwelling  distorts  the  near  surface  thermal  and  salinity  fields  by- 
transporting  colder,  more  saline  water  from  depth  up  to  the  surface,  which  typically 
results  in  density  fronts  along  the  coast.  During  periods  of  intense  upwelling,  the 
Davidson  Current  is  depressed  and  the  California  Current  may  be  displaced  as  much  as 
200km  away  from  shore  (Breaker,  1983). 

Bottom  topography  affects  current  flow,  upwelling,  and  tidal  flow  by  channeling 
the  flow  and  by  introducing  vorticity  into  the  flow.  The  major  bathymetric  feature  of 
the  Monterey  Bay  is  the  Monterey  Canyon.  The  canyon  has  a  northeast-southwest 
orientation  with  the  maximum  depth  ranging  from  100m  near  shore  to  3600m  at  the 
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continental  slope.  Within  the  study  area  the  canyon  axis  is  somewhat  aligned  with  the 
XE-SW  diagonal  of  the  study  area  with  depth  ranges  from  200m  at  northwest  and 
southeast  corners  to  2000m  at  the  southwest  corner  (Fig.  1.2).  Coastal  upwelling  has 
been  observed  to  be  greatly  intensified  over  the  heads  of  some  submarine  canyons  due 
to  the  channeling  effect  (Breaker,  19S3).  Additionally,  due  to  the  large  topographic 
variation  in  the  bathymetry  current  flow  perpendicular  to  the  axis  of  the  canyon  will 
tend  to  be  deflected  cyclonically  due  to  an  increase  in  positive  relative  vorticity  derived 
from  the  increase  (stretching)  of  the  water  column  as  the  flow  crosses  over  the  canyon 
axis.  Flow  parallel  to  the  canyon  axis  will  tend  to  be  channelled  along  the  axis. 
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R/V  RCRN1R  Cruise  Track 

(3  OCT  -  5  OCT  1981) 


Fig.  1.1  U.S.S.  DOLPHIN  and  R/V  ACANIA  operation  area. 
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122  24  W  121 °54  W 


Fi"  1  2  Bathvmetrv  of  the  Monterev  Bav  in  the  vicinity  of  the  study  area 
*  '  '  (contours  in  fathoms). 
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II.  THEORY 


Water  column  stability  is  an  expression  of  the  likelihood  for  the  occurance  of 
vertical  motion  of  water  parcels  within  the  water  column.  It  is  a  function  of  the 
density  variation  with  depth  (static  stability)  and  the  vertical  variation  in  velocity 
(shear).  Together,  the  magnitude  of  these  two  components  provide  an  indication  of  the 
dynamic  stability  of  the  water  column.  Irreversible  vertical  motion  (mixing)  of  water 
parcels  can  occur  when  there  is  an  imbalance  between  the  horizontal  velocity  shear  and 
the  thcrmohaline  stratification.  The  result  of  mixing  is  to  reduce  the  vertical  variation 
of  density  thereby  making  the  density  of  the  water  column  more  homogeneous.  This 
results  in  an  increase  in  the  center  of  gravity  of  the  water  column  and  a  consequent 
increase  in  the  potential  energy  of  the  system.  The  increase  in  potential  energy  is 
derived  from  the  kinetic  energy  of  the  vertical  motion  of  the  parcels  which  in  turn  is 
derived  from  the  kinetic  energy  of  the  force  causing  the  imbalance  in  the  density  or 
velocity  profiles  (Pond  and  Pickard,  1983,  p.60).  Definition  of  the  parameters  that 
define  stability  as  applied  to  this  study  are  discussed  below. 

A.  STATIC  STABILITY 

Static  stability  is  strictly  a  function  of  the  vertical  variation  in  the  density  field 
and  as  indicated  by  the  word  'static'  implies  no  initial  motion  in  the  system.  Stability 
is  indicated  by  the  variation  in  density  with  depth.  A  stable  density  profile  is  one  in 
which  density  increases  with  depth.  Consequently  gravitational  forces  will  not  impart 
any  motion  to  the  water  parcels.  However,  if  the  density  profile  decreases  with  depth 
then  heavier  water  overlies  lighter  water  and  the  result  is  for  the  lighter  water  to  rise 
while  the  heavier  water  sinks  (vertical  motion)  in  order  to  reestablish  a  gravitational 
balance. 

Density  in  the  ocean  is  a  function  of  temperature,  salinity,  and  pressure  in  order 
of  decreasing  importance.  However  a  more  useful  parameter  is  sigma-t  which  is 
defined  as: 

<7t  =  (p(t,s,0)  -  1000)  kg/m  3  (eqn  2.1) 
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1'his  parameter  is  useful  in  that  it  allows  for  a  better  estimate  of  what  the  density 
difference  between  two  water  parcels  would  be  if  they  were  at  the  same  level.  Sigma-t 
values  were  calculated  from  the  temperature  and  salinity  profiles  obtained  from  the 
CTD  data.  The  quantitative  evaluation  of  water  column  static  stability  is  given  by  the 
stability  parameter  "E"  which  is  expressed  as  follows: 


E  =  1  p  5<rt  5z 


(eqn  2.2) 


The  5(T(  Oz  term  represents  the  vertical  gradient  of  density. 

The  buoyancy  frequency  defines  the  natural  or  resonant  frequency  that  a  parcel 
would  oscillate  at  if  displaced  vertically  from  its  stable  position  within  a  statically 
stable  fluid.  It  is  quantified  as  the  square  root  of  the  product  of  the  static  stability 
parameter  "E"  and  the  gravitational  constant  "g"  and  is  symbolized  by  "N". 

X2  =  g  x  E  (eqn  2.3) 

The  buoyancy  frequency  is  also  the  maximum  frequency  of  internal  waves  in  water  of 
stability  "E"  (Pond  and  Pickard.  1983,  p.30). 

B.  DYNAMIC  STABILITY 

Dynamic  stability  is  an  indicator  of  the  effect  that  the  vertical  variation  in 
velocity  (shear)  will  have  on  the  generation  of  vertical  motion  turbulence  within  a 
water  column  of  stability  "E".  Therefore  it  is  a  function  of  both  the  static  stability  and 
the  magnitude  of  the  shear.  A  quantitative  expression  for  dynamic  stability  is  the 
"gradient  "  Richardson  number  which  is  the  ratio  of  the  square  of  the  buoyancy 
frequency  to  the  square  of  the  velocity  shear: 


Ri  =  (N2)  (5u.  5z)2 


(eqn  2.4) 


A  critical  value  of  0.25  has  been  established  from  turbulence  theory  as  the  cutoff 
between  dynamic  stability  and  instability  (Pond  and  Pickard,  19S3.  p.60).  Values  of  Ri 
>  14  indicate  that  the  vertical  gradient  of  velocity  is  insufficient  to  generate  sustained 
turbulence  and  therefore  dynamic  stability  is  is  sustained.  The  converse  is  true  for 
values  below  the  critical  value.  This  ratio  may  also  be  viewed  as  the  ratio  of  the 
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potential  energy  of  the  column  expressed  by  the  density  field  and  the  kinetic  energy 
expressed  by  the  velocity  field. 

C.  DOUBLE  DIFFUSION 

Instability  in  an  otherwise  statically  stable  environment  may  also  occur  from 
double  diffusion  which  is  the  result  of  the  difference  in  the  molecular  diffusivities  of 
temperature  and  salt  within  a  fluid  causing  an  imbalance  in  the  density  field.  The 
diffusivitv  of  temperature  is  approximately  100  times  greater  than  that  of  salt;  thus 
double  diffusion  may  occur  when  two  masses  of  approximately  the  same  density,  but 
different  temperature  and  salinity  are  in  contact  within  a  statically  stable  environment. 
The  difference  in  the  speed  of  molecular  diffusion  for  temperature  and  salt  results  in  a 
relatively  rapid  exchange  of  heat  without  a  compensating  exchange  of  salinity  thereby 
resulting  in  a  change  in  the  densities  of  the  two  water  masses.  There  are  three  possible 
combinations  of  temperature  and  salt  between  the  two  water  masses  which  result  in 
three  entirely  different  conditions:  ie.  double  diffusive  instability,  layering,  and  no 
motion. 

1.  Case  I:  Double  Diffusive  Instability 

Double  diffusive  instability  may  occur  if  there  is  warmer,  saltier  water  over 
cooler,  fresher  water  (negative  temperature  and  salinity  gradients).  The  upper  layer 
cools  while  the  lower  layer  warms  however,  the  salinity  exchange  is  not  rapid  enough 
to  compensate  for  the  increasing  density  of  the  upper  layer  and  the  decreasing  density 
of  the  lower  layer.  Thus,  a  negative  density  gradient  is  formed,  and  the  boundary 
between  the  two  layers  becomes  unstable.  Consequently,  the  lower  layer  tries  to  rise 
while  the  upper  layer  tries  to  sink  thereby  resulting  in  mixing.  Experiments  have 
shown  this  mixing  to  occur  as  thin  columns/filaments  of  each  layer  infiltrating  the 
adjacent  layer,  this  type  of  mixing  is  described  as  "salt  fingering"  (Pond  and  Pickard, 
19S3,  p. 3 1 ) 

2.  Case  II:  Layering 

An  intensification  of  the  stable  density  gradient  between  layers  occurs  when 
cooler,  fresher  water  is  over  warmer,  saltier  water  (positive  temperature  and  salinity 
gradients).  As  the  upper  layer  warms  the  density  decreases  while  the  converse  is 
happening  in  the  lower  layer.  Again  the  salinity  exchange  is  not  rapid  enough  to 
compensate  for  the  density  changes.  When  the  temperature  difference  is  large  the 
warming  of  the  upper  layer  by  the  lower  layer  will  cause  it  to  rise  due  to  its  reduced 
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density  .  The  lower  layer  will  attempt  to  sink  due  to  cooling  and  increased  density. 
This  will  lead  to  a  breakdown  in  the  interface  between  the  layers  thereby  resulting  in 
mixing. 

3.  Case  III:  No  Motion 

In  the  case  of  warmer,  fresher  water  over  cooler,  saltier  water  (negative 
temperature  gradient,  positive  salinity  gradient)  the  density  difference  is  reduced  as  the 
upper  layer  cannot  get  denser  than  the  layer  beneath  it.  Therefore  no  motion  occurs. 
In  order  for  double  diffusion  to  result  in  instability,  the  gradients  of  temperature  and 
salt  must  be  the  same  (Pond  and  Pickard.  19S3,p.32). 


16 


III.  DATA  ACQUISITION  AND  PROCESSING 


A.  DATA  ACQUISITION 

Temperature,  conductivity,  and  velocity  data  analyzed  in  this  study  were  acquired 
from  3  to  5  October,  1984  by  the  L'.S.S.  DOLPHIN  and  R/V  ACANIA,  which 
sampled  a  square-shaped  region  centered  at  36°  44.5'N,  122°  03.0'W,  approximately  22 
km  west  of  Moss  Landing.  California  (Fig.  1.1).  The  R/V  ACANIA  made  seven 
circuits  of  the  edge  of  the  region  while  U.S.S.  DOLPHIN  sampled  the  interior  of  the 
region.  The  vessels  were  equiped  as  follows. 

I.  U.S.S.  DOLPHIN  Instrumentation 

The  U.S.S.  DOLPHIN  was  fitted  with  a  5m  high  tripod  structure  near  the 
bow  of  the  submarine.  This  structure  was  used  to  support  the  temperature, 
conductivity,  and  depth  sensors.  The  CTD  sensors  consisted  of  Seabird  Electronics 
thermistor  and  conductivity  cell,  and  a  Viatran  model  304  strain  gauge  pressure  sensor. 
The  data  from  these  sensors  were  digitized  and  stored  on  9-track  magnetic  tape. 

Velocity  profiles  relative  to  the  submarine  were  measured  by  T.P.  Stanton  with 
an  R-D  Instruments  acoustic  doppler  profiler  mounted  at  deck  level,  2m  forward  of  the 
base  of  the  instrument  tripod.  The  acoustic  profiler  had  an  upward-looking,  4-beam 
Janus  configuration.  The  beam  angles  were  rotated  45°  in  azimuth  from  the  forward 
and  cross  axes  of  the  submarine  in  order  to  minimize  acoustic  reflection  from  the 
instrument  tripod.  Despite  this  rotation,  reflections  from  the  tripod  were  received  by 
the  aft  transducer  sidelobes  causing  the  second  and  third  velocity  bins  to  have 
unreliable  data.  The  acoustic  beams  were  inclined  30°  from  the  vertical  axis.  Acoustic 
pulses  were  transmitted  at  1.2  MHz  from  all  four  beams  at  a  repetition  rate  of  5 
pings  sec.  The  received  signal  was  range-gated  into  30  sections  corresponding  to  lm 
length  vertical  bins,  providing  profiles  of  doppler  velocities  in  the  30m  above  the  hull  of 
the  submarine.  The  radial  doppler  shifts  were  resolved  into  orthogonal  velocity 
components  by  differencing  the  doppler  shifts  of  each  bin  in  opposing  beams  using  the 
relationship: 

Vl3  =  ((f\  -  fj)  4fg)  x  CQ  sin(0)  (eqn  3.1) 
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V24  =  <(f2  ’  f4>  Jf0>  X  C0  Sln(G) 


(eqn  3.2) 


V..:  velocity  component  resolved  from  the  "i"th  and  "j "th  beams 
fQ:  transmitted  frequency 

f:  received  doppler  shifted  frequency  from  the  "i"th  beam 
CQ:  velocity  of  sound  at  the  transducer  face 
0:  beam  inclination  from  the  vertical 

These  two  orthogonal  velocity  measurements  were  then  rotated  by  45°  to  give  the 
forward  and  cross  velocities.  Velocity  components  were  averaged  for  10  seconds  to 
reduce  the  high  noise  level  inherent  in  pulsed  incoherent  acoustic  doppler  systems 
which  principally  arises  from  trying  to  resolve  small  doppler  frequency  shifts  within 
short  range  intervals.  The  platform  referenced  velocity  components  were  recorded  and 
displayed  on  a  Hewlett-Packard  216  computer. 

Other  instrumentation  onboard  R  V  DOLPHIN,  though  not  used  in  this 
study,  included  microstructure  measuring  equipment  used  by  T.R.  Osborn  and  R.G. 
Lueck.  These  consisted  of  3  fast  response  thermistors,  2  airfoil  shear  probes  and  3-axes 
accelerator.  The  airfoil  probes  were  used  to  measure  the  time  derivative  of  orthogonal 
components  of  horizontal  velocity  (du  dt,  dv/dt).  The  time  derivative  is  then  converted 
into  a  measurement  of  shear.  The  scales  of  the  shear  measured  by  this  method  is  on 
the  order  of  0.01m  to  0.5m.  Use  of  these  instruments  is  described  in  detail  by  Osborn 
and  Lueck  (1985).  The  fast  response  thermistors  were  used  to  measure  fine-structure 
temperature  gradients.  The  analysis  of  these  measurements  of  mixing  activity  is  being 
performed  by  T.R.  Osborn  and  E.C.  Itsweire. 

Additionally,  there  were  two  acoustic  transducers  which  transmitted  an  80 
KHz  pulse  in  the  forward  and  vertical  directions  with  the  intent  of  measuring  the 
scattering  cross-section  of  the  backscatterers  in  the  volume  ahead  and  above  the 
submarine.  These  were  used  by  D.  Farmer  and  C.  Crawford  to  investigate  the  nature 
of  acoustic  back  and  side  scatter  at  this  frequency. 

2.  R/V  ACANIA  Instrumentation 

R.'V  ACANIA  utilized  several  types  of  instruments  to  measure  the 
temperature  and  conductivity  of  the  upper  ocean  between  the  surface  and  in  excess  of 
115m  depth.  Ocean  skin  temperature  was  measured  with  a  Rosemont  platinum 
resistance  thermometer  (5.0  mdeg  absolute  calibration)  suspended  from  a  boom  over 
the  side  of  the  ship  such  that  the  thermometer  remained  within  the  upper  few 
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centimeters  of  the  surface.  Near  surface  (2m)  temperatures  and  conductivities  were 
measured  with  a  pump-through  Seabird  thermistor  and  conductivity  cell  located  in  the 
ship's  seachest.  Continuous  measurement  of  temperature  and  conductivity  between  5 
and  115  meters  was  made  by  towvoing  a  Neil  Brown  Instrument  Systems  Mk  111  CTD. 
The  CTD  consisted  of  a  combination  thermistor  and  platinum  resistance  thermometer 
(0.5  mdeg  resolution,  5.0  mdeg  absolute  calibration),  platinum  electrode  conductivity 
cell  (0.001  mmho  resolution,  0.01  nominal  accuracy),  and  an  electrical  strain  gauge 
pressure  sensor  (10  cm  resolution).  During  the  towyo  the  CTD  package  was  raised  and 
lowered  between  specified  depth  limits  while  maintaining  a  constant  ship  speed  of 
approximately  5  knots.  By  this  method  a  continuous  series  of  near  vertical 
temperature  and  conductivity  profiles  were  obtained.  A  spatial  resolution  for  the 
vertical  section  was  achieved  on  the  order  of  0.5km  by  0.5m  based  on  the  distance 
covered  in  the  time  to  complete  an  upcast/downcast  and  the  binning  of  temperature 
and  conductivity  every  0.5m  of  depth.  This  sampling  resolution  was  used  to  determine 
the  T,  S,  <?t.  and  N  structure  of  the  upper  100m  of  the  water  column  on  a  scale 
comparable  with  the  shipboard  ADVP  measurements.  Thus  a  valid  comparison 
around  the  perimeter  of  the  study  area  could  be  made  between  the  static  stability  as 
determined  from  the  CTD  data  and  the  shear  stability  determined  from  the  ADVP 
data.  Location  (latitude/  longitude)  information  for  each  CTD  profile  was  obtained 
from  an  Internav  LC408  LORAN  C  receiver  (25m  resolution,  100m  accuracy). 

Velocity  profile  measurements  below  the  ship  were  made  with  an  Ametek- 
Straza  acoustic  doppler  profiler  model  DCP4015  (ADVP).  This  profiler  has  a  4-beam 
Janus  configuration  transducer  aligned  with  the  forward  and  cross  axes  of  the  ship 
(Fig.  3.1).  Acoustic  pulses  were  transmitted  at  300KHz  every  0.6  seconds  (1.67 
pings  sec)  with  a  5  msec  pulse  duration.  The  received  signal  was  range-gated  into  32 
bins  giving  a  3m  vertical  bin  size;  the  shallowest  range  bin  was  centered  at  7.5m  depth. 
The  forward  and  cross  velocity  components  for  each  depth  level  were  computed  by 
differencing  doppler  shifts  from  complementary  beams  as  in  the  U.S.S.  DOLPHIN 
(section  A.).  Forward  and  cross  velocity  estimates  from  each  ping  were  rotated  into  a 
North-South  coordinate  system  using  heading  information  from  the  ship's  Mk9  Sperry 
gyro.  The  resolved  velocity  components  were  then  averaged  over  2  minutes  in  order  to 
reduce  the  doppler  frequency  noise  level  as  well  as  noise  due  to  the  pitch  and  roll  of 
the  platform.  As  indicated  by  Kosro  (1985),  when  using  comparable  averaging 
intervals  and  moderate  sea  conditions,  pitch  and  roll  affects  the  vertical  velocity 
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component  the  most  with  velocity  overestimates  on  the  order  of  3  -  4  cm  sec. 
However,  the  horizontal  velocity  components  are  affected  on  the  order  of  less  than  1 
cm  sec:  thus  the  horizontal  velocity  estimates  may  be  calculated  with  minimal  error 
when  pitch  and  roll  corrections  are  ignored,  particularly  for  the  low  sea  states 
encountered  during  these  measurements. 

The  ADVP  measures  velocity  profiles  relative  to  the  vessel.  Absolute  current 
profiles  may  be  derived  by  measuring  and  removing  the  ground  referenced  ship's 
velocity.  In  this  study  only  the  shear  profiles  will  be  used  as  high  resolution  profiles 
were  needed  for  the  construction  of  Richardson  Number  profiles  and  the  ship's  velocity 
cannot  be  determined  to  comparable  accuracy  with  only  two  minutes  averages: 
typically  15  minutes  are  needed  to  resolve  the  ship's  velocity  to  10  cm  sec  using 
LORAX  C. 

B.  DATA  PROCESSING 

The  stability  of  the  water  column  and  the  processes  affecting  it  were  examined  by 
analyzing  the  temperature,  salinity,  sigma-t,  and  velocity  fields.  Stability  is  a  function 
of  the  vertical  gradients  of  density  and  velocity,  but  since  the  water  mass  structure  is 
not  readily  apparent  in  the  density  field  alone,  the  temperature  and  salinity  fields  were 
also  analyzed.  Analysis  was  accomplished  by  the  construction  of  graphical  displays  of 
contoured  vertical  sections,  profile  plots,  and  horizontal  vector  plots.  The  temperature, 
salinity,  sigma-t,  and  static  stability  fields  were  constructed  from  the  R  V  ACANIA 
CTD  data,  which  provided  information  from  5.0m  to  115.0m  depths.  These  fields  are 
displayed  as  both  contoured  vertical  sections  and  profile  plots.  U.S.S.  DOLPHIN 
CTD  data  was  not  used  in  the  construction  of  these  fields  as  information  was  limited 
to  the  cruise  depth  of  the  submarine  and  thus  had  limited  vertical  range.  Velocity 
fields  were  constructed  from  both  the  R  V  ACANIA  and  L’.S.S.  DOLPHIN  ADVP 
data  sets  and  were  displayed  as  profile  plots  and  horizontal  vector  plots.  Richardson 
number  fields  were  constructed  by  combining  the  sigma-t  and  velocity  fields  from  the 
R  V  ACANIA  CTD  and  ADVP  data  sets  and  are  displayed  as  profile  plots. 

Contoured  vertical  sections  were  chosen  for  the  temperature,  salinity,  and  sigma-t 
fields  because  they  provided  a  ready  depiction  of  the  spatial  variability  of  a  parameter 
within  a  particular  section  as  well  as  an  indication  of  the  temporal  variability  through 
successive  vertical  sections.  The  vertical  sections  were  constructed  from  2.0m  vertical 
averages  of  the  R  V  ACANIA  CTD  data.  Vertical  averaging  was  necessary  in  order  to 
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reduce  the  noise  level  in  the  contouring  procedure.  Individual  vertical  sections  may  be 
considered  to  be  approximate  "snapshots"  of  the  variability  expressed  in  the 
temperature,  salinity,  and  sigma-t  fields  as  they  represent  data  collected  within 
approximately  one  hour  along  a  Skm  horizontal  section. 

Averaged  and  staggered  profile  plots  were  constructed  for  the  parameters  of 
temperature,  salinity,  sigma-t,  velocity,  static  stability,  and  Richardson  number.  The 
plots  are  useful  in  that  they  provide  a  more  detailed  depiction  of  the  vertical  variation 
of  a  parameter  than  does  the  contoured  vertical  section.  Thus,  small  scale  vertical 
variability  associated  with  some  processes  such  as  internal  waves  and  double  diffusion 
are  more  readily  apparent.  Also,  in  combination  with  the  contoured  vertical  sections,  a 
better  understanding  of  the  temporal  and  spatial  variability  of  the  fields  of  interest  may 
be  made.  The  averaged  profiles  of  temperature,  salinity,  and  sigma-t  were  constructed 
by  side-averaging  the  vertical  profiles  of  these  parameters  at  2.0m  vertical  resolution. 
The  standard  deviation  about  these  profiles  was  also  calculated  and  depicted  by  dashed 
profiles  on  either  side  of  the  averaged  profile.  Included  with  the  side-averaged  profiles 
is  a  composite  profile  which  represents  the  mean  of  the  side-averaged  profiles. 
Averaged  profiles  of  static  stability  were  constructed  from  3.0m  vertical  resolution  side- 
averaged  profiles  to  correspond  with  the  resolution  of  the  ADVP  data  used  in  the 
calculation  of  the  Richardson  number;  standard  deviation  was  not  calculated.  Critical 
value  limits  are  expressed  on  the  E  and  Ri  plots  by  a  vertical  dotted  line.  Staggered 
profiles  were  constructed  for  temperature,  salinity,  sigma-t,  and  static  stability  from  the 
0.5m  resolution  CTD  data. 

Horizontal  vector  plots  of  velocity  were  constructed  for  the  25m  and  75m  depth 
levels  from  the  R,  V  ACANIA  ADVP  data  for  the  second  through  seventh  transits  of 
the  study  area.  The  velocity  difference  between  these  two  levels  (shear)  was  also 
constructed.  Together  these  plots  served  to  provide  information  on  the  tidal  and  large- 
scale  current  flow.  The  25m  and  75m  depth  level  was  chosen  because  they  defined  the 
upper  and  lower  depth  limits  of  the  section  of  the  water  column  exhibiting  the  greatest 
activity  in  terms  of  temperature  and  salinity  variability.  This  permitted  a  comparison 
between  the  dynamic  (velocity)  forcing  and  the  static  forcing  on  stability. 


DOPPLER  PROFILER  BEAM  CONFIGURATION 


Fig.  3.1  Illustration  of  4-beam  Janus  configured  ADVP. 
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STARBOARD 


IV.  OBSERVATIONS 


A.  GENERAL  DESCRIPTION 

A  general  depiction  of  the  spatial  variability  present  in  the  analyzed  fields  of 
temperature,  salinity  and  sigma-t  are  presented  in  the  side-averaged  contoured  vertical 
sections  (Figs.  4.1  -  4.3).  It  is  readily  apparent  that  there  is  considerably  more 
variability  in  the  temperature  and  salinity  vertical  sections  than  in  the  density  sigma-t 
vertical  section.  This  is  due  to  a  tendency  for  the  variations  in  temperature  and  salinity 
to  compensate  each  other  resulting  in  a  more  consistent  density  field.  This  is  evidenced 
particularly  by  the  compensation  exhibited  in  the  eastern  side  panel  comparing  the 
strong  thermal  gradient  between  35m  and  45m  with  the  salinity  minimum  at  the  same 
depth  levels.  Overall,  the  density  field  strongly  reflects  the  temperature  field;  however, 
where  there  are  strong  salinity  gradients,  salinity  variations  are  reflected  in  the  density 
field. 

Mean  profiles  of  temperature,  salinity,  and  sigma-t  for  each  side  and  transit  of 
the  study  area  depict  the  temporal  variability  in  the  vertical  variation  of  these 
parameters  (Figs.  4.4  -  4.15).  Of  note  are  the  consistently  strong  thermal  gradients 
between  15m  and  35m  depicted  in  the  profiles  for  the  east  and  south  sides  (Figs.  4.5 
and  4.6).  The  north  and  west  sides  exhibit  considerably  more  temporal  variability  in 
their  temperature  profiles  (Figs.  4.4  and  4.7).  This  variability  is  oscillatory  such  that 
the  composite  profile  (average  of  all  seven  circuits  of  a  side)  has  nearly  constant  slope. 
The  salinity  profiles  (Figs.  4.8  -  4.11)  suggest  a  migration  of  parcels  of  less  saline  water 
between  25m  and  65m  through  the  study  area.  Sigma-t  profiles  for  the  east  and  south 
sides  (Figs.  4.13  and  4.14)  reflect  the  temperature  profiles  especially  with  regard  to  the 
strong  thermal  gradients  as  was  indicated  in  the  contoured  vertical  sections.  The 
sigma-t  profiles  for  the  north  and  west  sides  (Figs.  4.12  and  4.15)  however  do  not  tend 
to  reflect  the  vertical  variability  evident  in  the  temperature  profiles,  and  instead  tend  to 
monotonically  increase  with  depth.  This  may  be  attributed  to  compensation  by  the 
salinity  field. 

B.  STATIC  AND  DYNAMIC  FORCING  MECHANISMS 

As  discussed  in  Chapter  I,  there  are  several  forcing  mechanisms  which  can  induce 
variability  in  the  temperature,  salinity,  density,  and  velocity  fields  and  thus  ultimately 
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the  stability  of  the  water  column.  The  mechanisms  investigated  in  this  study  were 
upwelling,  bottom  topography,  tidal  flow,  current  systems,  and  internal  gravity  waves. 
Analysis  of  the  temperature,  salinity,  sigma-t.  and  velocity  fields  in  conjunction  with 
the  atmospheric  and  tidal  conditions  present  during  the  study  indicated  that  all  of  these 
mechanisms  were  present  with  the  exception  of  upwelling.  Evidence  for  these 
mechanisms  is  discussed  below. 

1.  Upwelling 

The  atmospheric  conditions  several  days  prior  to  and  during  the  study  period 
were  characterized  by  a  high  pressure  ridge  off  the  central  California  coast  with  a  weak 
pressure  gradient  over  the  coastal  areas.  The  resulting  low  level  winds  were  from  the 
north-northwest  at  approximately  5  -  8  knots.  This  weak  NNW  flow  was  modified  by 
the  diurnal  land-sea  breeze  fluctuations.  A  sea  breeze  regime  commenced  in  the  early 
afternoon  with  a  backing  and  increase  in  wind  speed  to  WNW  at  8  -  14  knots.  This 
regime  persisted  until  midnight  when  the  winds  veered  back  to  NNW  and  diminished  to 
2  -  4  knots  (weak  land  breeze  regime).  The  synoptic  scale  flow  resumed  in  the  early 
morning  just  prior  to  sunrise.  The  low  winds  and  generally  offshore  flow  resulted  in  an 
absence  of  the  summertime  pervasive  coastal  status  and  fog  .  Low  winds  also  provided 
for  calm  sea  conditions.  Due  to  the  weak  NNW  flow,  coastal  upwelling  was  severly 
diminished  and  thus  deemed  to  be  negligible  as  a  forcing  mechanism  during  the  study 
period. 

2.  Bottom  Topography 

The  bottom  topography  of  the  study  area,  as  mentioned  previously,  was 
dominated  by  the  steep  slopes  of  the  Monterey  Canyon  (Fig.  1.2)  The  canyon  axis 
roughly  corresponded  with  the  NE-SW  diagonal  of  the  study  area  with  slope  towards 
the  southwest,  but  is  deflected  somewhat  from  the  NE-SW  orientation  by  the  presence 
of  a  large  topographic  ridge  projecting  southward  into  the  study  area  from  near  the 
midpoint  of  the  north  side.  Evidence  of  the  effects  of  bathymetry  on  the  analyzed 
fields  is  demonstrated  in  the  contoured  vertical  sections  of  temperature  (Figs.  4.16  - 
4.22).  Particularly  evident  is  the  mass  of  colder  water  between  40m  and  95m  in  the 
eastern  side.  This  mass  is  expressed  in  the  vertical  section  by  an  upward  bulging  of  the 
11.5  and  12.0°C  isotherms  in  the  northern  half  of  the  eastern  side.  The  bathymetry  of 
the  eastern  side  is  nearly  constant  at  approximately  400  fathoms  (800m)  for  the 
northern  thrce-forths  of  the  side  then  rises  to  100  fathoms  (200m)  at  the  southeast 
corner.  Therefore,  this  colder  mass  would  be  indicative  of  a  trapped  pool  of  colder 
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water.  Analysis  of  the  12°C  isotherm  in  the  northern  side  sections  also  show  evidence 
of  bathymetric  channeling  as  expressed  by  an  downward  bulging  in  the  isotherm.  The 
location  of  the  bulging  corresponds  with  the  axis  of  the  canyon  and  would  seem  to 
represent  the  channeling  of  warmer  water  along  the  canyon  axis.  The  combination  of 
the  warmer  pool  in  the  eastern  portion  of  the  northside  with  the  cooler  pool  in  the 
northern  portion  of  the  eastern  side  indicates  that  warmer  water  is  being  advected 
along  the  canyon  axis  while  the  cooler  water  is  displaced  along  the  sides  of  the  canyon. 
Channeling  of  water  flow  by  the  canyon  may  also  affect  the  velocity  field,  especially 
with  regard  to  the  oscillatory  tidal  flow. 

3.  Tidal  Flow 

Tidal  flow  along  the  California  Coast  is  characterized  by  two  high  and  two 
low  tides  every  25  hours,  occurring  about  50  minutes  later  every  day.  The  heights  of 
the  two  high  tides  and  the  heights  of  the  two  low  tides  are  not  equal  in  magnitude  or  in 
time  between  them.  During  the  study  period  the  higher  high  tide  occurred  in  the  early 
evening  at  about  1900  while  the  lower  low  tide  occurred  shortly  after  midnight.  Actual 
tide  times  at  Moss  Landing  were  high  tide  at  ISIS,  3  Oct  and  0839  and  1924,  4  Oct; 
low  tide  occurred  at  0140  and  1403,  4  Oct  and  0240.  5  Oct. 

Tidal  effects  are  most  noticeable  on  the  temperature  and  velocity  fields.  This 
is  demonstrated  in  the  contoured  vertical  sections  of  temperature  (Figs.  4.16  -  4.22)  and 
the  horizontal  vector  plots  of  velocity  (Figs.  4.23  -  4.27).  With  regard  to  the 
temperature  field,  tidal  effects  are  especially  evident  in  the  vertical  sections  for  the 
northern  side  expressed  by  the  variability  in  the  12°  and  14°C  isotherms.  This  may  be 
due  to  the  near  parallel  orientation  of  the  section  with  the  predominantly  east-west 
fluctuations  in  the  near  surface  tidal  flow.  In  circuit  1  (Fig.  4.16)  the  12°  isotherm  is 
nearly  fiat  across  the  northern  section  at  a  depth  of  approximately  60m  during  a  time 
corresponding  to  a  period  of  slack  water  at  high  tide.  By  circuit  2  (Fig.  4.17),  the  12° 
isotherm  slopes  upward  to  the  east  and  a  patch  of  warmer  water  is  evident  in  the 
eastern  portion  of  the  section  between  60  and  80m;  this  is  during  a  period  of  ebb  flow. 
In  circuit  3  (Fig.  4.18),  a  period  of  flood  flow,  the  12°  isotherm  is  deeper  and  slopes  in 
the  opposite  direction  to  the  previous  circuit,  but  exhibits  a  bulge  in  the  isotherm  to 
representing  an  intrusion  of  warmer  water.  At  low  tide  (Fig.  4.20)  the  12°  isotherm  is 
again  nearly  fiat  as  in  circuit  1,  however  the  depth  is  nearly  15m  deeper  indicating  the 
ebb  flow  introduced  warmer  water  into  the  study  area  which  depressed  the  12' 
isotherm.  The  tilt  of  the  14°  isotherm  also  reflects  the  tidal  flow.  During  periods  of 
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ebb  flow  (Fig.  4.17  and  4.22).  the  isotherm  slopes  downward  to  the  east  whereas  the 
slope  is  reversed  for  flood  flow  (Fig.  4.18  and  4.21).  The  14°  isotherm  is  approximately 
horizontal  during  periods  of  slack  water  (f  ig.  4.16  and  4.20).  Tidal  flow  may  also 
distort  the  temperature  and  salinity  fields  through  the  advection  of  water  parcels  such 
that  the  advected  parcels  have  different  temperature  and  salinity  characteristics  than 
the  host  water.  Evidence  of  tidal  induced  advection  was  discussed  above  with  regard  to 
the  influence  of  bottom  topography. 

The  horizontal  velocity  plots  (Figs.  4.23  -  4.27)  exhibit  fluctuations  in 
direction  and  magnitude  corresponding  to  the  tidal  period.  This  is  particularly 
evidenced  by  the  eastern  and  western  sides.  The  western  side  shows  a  consistent 
westerly  flow  .  however  the  magnitude  varies  in  direct  relationship  to  the  east-west  flow 
reversals  observed  on  the  eastern  side.  The  period  between  the  flow  reversals  is  on  the 
order  of  13  hours  which  corresponds  well  with  the  tidal  period.  The  northern  and 
southern  sides,  on  the  other  hand,  exhibit  consistent  southerly  flow  through  the  study 
period:  this  could  therefore  be  indicative  of  a  general  southerly  flow  (California 
Current)  with  a  tidally  induced  fluctuating  east-west  velocity  component.  Bottom 
topography  induced  channeling  results  in  the  flow  fluctuations  being  greater  along  the 
eastern  side. 

4.  Currents 

Current  systems  of  interest  to  this  study  included  the  California  and  Davidson 
Currents.  The  Davidson  Current  is  usually  dominant  during  periods  of  reduced 
upwelling  as  occurred  during  this  study.  However,  the  horizontal  vector  velocity  plots 
(Figs.  4.23  -  4.27)  along  the  northern  and  southern  sides  of  the  study  area  do  not  give 
any  indication  of  the  presence  of  the  northward  flowing  Davidson  Current  within  the 
study  area,  but  rather  depict  a  consistent  southward  flow  characteristic  of  the 
California  Current.  As  mentioned  above  the  flow'  along  the  eastern  and  western  sides 
appears  to  be  modified  by  tidal  flow.  Additionally,  the  current  velocities  are  greater  at 
25m  than  at  75m  resulting  in  a  negative  velocity  gradient  and  providing  for  shear  as  a 
dynamical  forcing  on  stability.  Average  velocity  differences  between  the  25m  and  75m 
depth  levels  was  on  the  order  of  5  cm  sec. 

Currents  also  affect  the  temperature  and  salinity  fields  through  advection  of 
water  parcels.  This  appears  evident  in  the  contoured  vertical  sections  of  salinity  (Figs. 
4.2S  -  4.34)  where  a  lens  shaped  minimum  in  salinity  is  described  by  the  33.3  ppt 
isohaline  at  between  25m  and  55m  depth.  Advection  of  this  less-saline  lens  by  the 
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large-scale  current  flow  may  be  observed  in  the  vertical  sections  for  the  north  and  east 
sides  where  the  lens,  which  initially  spans  the  northside  at  approximately  45m,  is 
observed  in  successive  sections  to  reduce  in  size  and  migrate  to  the  east  and  south. 
This  southward  migration  is  also  evident  in  the  vertical  sections  for  the  south  and 
westsides  where  smaller  parcels  of  less-saline  water  are  observed  to  migrate  across  the 
side  during  the  course  of  the  study.  An  estimate  of  the  speed  of  migration  is 
approximately  0.5  km.  hr  (14cm/s).  This  is  only  a  crude  estimate  considering  the  other 
processes  which  may  distort  the  parcels.  The  overall  southward  migration  apparent  in 
these  vertical  sections  is  consistent  with  the  general  southerly  flow  observed  in  the 
velocity  plots  and  may  be  attributed  to  the  California  Current  system. 

5.  Internal  Gravity  Waves 

Internal  gravity  waves  affect  the  temperature,  salinity,  and  density  fields 
through  the  intensification  and  weakening  of  gradients  caused  by  wrave  induced 
displacements  of  the  density  interface.  The  movement  of  the  internal  wave  intensifies 
temperature,  salinity,  and  density  gradients  in  the  direction  of  movement  as  the  density 
interfaces  oscillates  vertically.  A  concurrent  weakening  of  gradients  occurs  behind  the 
movement  of  the  density  interface.  The  frequency  of  internal  waves  is  bounded  at  the 
upper  end  by  the  buoyancy  frequency.  A  mean  buoyancy  frequency  of  approximately 
0.6  eph  was  calculated  for  the  study  area  though  this  represents  a  low  value  for  the 
upper  bound  of  the  buoyancy  frequency  as  individual  profiles  values  were  higher.  This 
frequency  corresponds  to  a  period  of  1.7  hours  indicating  that  the  shorter  period 
internal  waves  should  be  adequately  represented  by  these  diagrams.  The  lower  bound 
on  the  frequency  for  internal  waves  is  the  inertial  frequency.  The  inertial  frequency  for 
the  latitude  of  the  study  area  (36.75°N)  was  calculated  to  be  0.0498  cph.  The 
corresponding  inertial  is  20.06  hours. 

While  wave-like  features  appear  evident  in  the  contoured  vertical  sections, 
they  are  especially  evident  in  the  staggered  profile  plots.  This  is  particularly  noticeable 
in  the  staggered  profile  plots  for  the  eastern  and  western  sides  (Figs.  4.35  -  4.37).  In 
the  staggered  temperature  profiles  for  the  eastside,  the  area  of  strongest  thermal 
gradient  (between  15m  and  35m  depth)  is  observed  to  oscillate  in  a  wave-like  pattern 
with  a  wavelength  of  approximately  7km.  A  second  wave  pattern  is  evident  between 
65m  and  85m  and  is  in  phase  with  the  upper  wave  pattern.  Similar  wave  patterns  are 
exhibited  in  the  staggered  profiles  of  salinity  and  sigma-t.  The  western  side  staggered 
temperature  profiles  exhibit  two  sets  of  waves  at  between  10m  and  30m  for  the  upper 


27 


wave  and  between  75m  and  95m  for  the  lower  wave  pattern.  These  two  wave  patterns 
are  ISO0  out  of  phase  which  results  in  a  stronger  alternating  intensification  and 
weakening  of  density  gradients  between  them.  As  in  the  staggered  profiles  for  the 
eastern  side,  the  western  side  salinity  and  sigma-t  profiles  exhibit  similar  wave  patterns 
to  that  displayed  in  the  temperature  profiles.  This  similarity  in  wave  patterns  between 
all  three  parameters  is  due  to  the  movement  of  internal  gravity  waves  along  density 
interfaces  resulting  in  an  oscillatory  or  wave  motion  expressed  in  the  density  profiles. 
Since  the  density  is  not  changed  by  the  motion  of  the  internal  wave,  the  temperature 
and  salinity  profiles  should  exhibit  similar  wave  patterns. 

C.  STATIC  AND  DYNAMIC  STABILITY 

The  variability  expressed  in  the  displays  of  temperature  and  salinity  were  for  the 
most  part  compensating  such  that  the  displays  of  sigma-t  exhibited  a  more  uniform 
field.  This  is  exemplified  in  the  contoured  vertical  section  of  sigma-t  (Figs.  4.38  -  4.44). 
Consequently,  the  static  stability  field  could  be  expected  to  be  stable  overall  as 
exhibited  by  displays  of  the  static  stability  parameter  (E)  and  the  dynamic  stability 
parameter  (Ri)  or  Richardson  number.  Comparison  of  the  averaged  profiles  of  the 
static  stability  parameter  (E)  (Figs  4.45  -  4.4S)  and  the  Richardson  number  (Ri)  (Figs. 
4.49  -  4.52)  provide  an  indication  of  the  relative  importance  between  the  vertical 
density  gradient  and  the  velocity  gradient  on  the  stability  of  the  water  column.  The  E 
and  Ri  average  profiles  are  calculated  from  3m  vertically  averaged  density  differences 
to  coincide  with  the  3m  vertical  resolution  of  the  velocity  data.  Consequently,  because 
gradients  of  density  and  shear  can  occur  on  scales  smaller  than  this  3m  differencing 
interval,  the  Ri  profile  values  represent  an  upper  bound.  Many  of  the  maximums  in 
the  Ri  profiles  correspond  with  the  maximums  in  the  E  profiles  indicating  that  the 
strength  of  the  density  gradient  is  driving  the  dynamic  stability.  Maximums  in  the  Ri 
profiles  which  do  not  correspond  to  maximums  in  the  E  profiles  indicate  that  there  is  a 
relative  minimum  in  the  vertical  gradient  of  velocity  (shear)  in  that  portion  of  the  water 
column  as  well  as  a  weak  density  gradient.  Minimums  in  the  Ri  profiles  indicate  that 
the  relative  magnitude  of  the  shear  is  greater  than  the  density  gradient.  Values  of  Ri 
below  the  critical  value  indicate  that  the  shear  is  greater  than  twice  the  density 
gradient.  Of  note  are  the  below  critical  Ri  values  observed  near  the  top  of  most 
profiles.  This  is  a  manifestation  of  the  acoustic  profiler  caused  by  sidelobe  interference 
in  the  upper  bins  of  the  shipboard  ADVP. 
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The  averaged  profiles  of  the  static  stability  parameter  indicate  stability  through 
the  vertical  extent  of  the  water  column.  However,  the  staggered  profiles  of  0  indicate 
thin  patches  ( <  2.0m  thick)  of  static  instability.  This  is  exemplified  in  the  staggered 
profiles  for  the  east  (Fig.  4.53)  and  west  sides  (Fig.  4.54).  As  was  noted  in  the 
discussion  of  internal  waves,  wave  like  patterns  are  also  exhibited  in  the  staggered  E 
profiles.  The  thin  unstable  patches  are  indicative  of  double  diffusive  processes  acting 
upon  the  disturbances  in  the  density  field  caused  by  the  static  and  dynamic  forcing 
mechanisms  of  bottom  topography,  tidal  and  current  flow,  and  internal  waves. 
Preferred  locations  of  instability  and  double  diffusive  processes  are  in  regions  of  strong, 
but  opposing  temperature  and  salinity  gradients.  However,  despite  the  extreme 
variability  observed  in  the  temperature  and  salinity  fields,  the  density  fields  were  for  the 
most  part  consistent  with  monotonically  increasing  gradients,  and  consequently  the 
water  column  is  predominantly  stable. 
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Fig.  4.1  Vertical  section  of  side-averaged  temperature. 
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Fig.  4.2  Vertical  section  of  side-averaged  salinity. 
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Fig.  4.3  Vertical  section  of  side-averaged  sigma-t. 
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Fig.  4.5  Averaged  profiles  of  temperature  (eastside). 
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Fig.  4.6  Averaged  profiles  of  temperature  (southside). 
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Fig.  4.7  Averaged  profiles  of  temperature  (westside). 
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Fig.  4.8  Averaged  profiles  of  salinity  (northside). 
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Fig.  4.9  Averaged  profiles  of  salinity  (eastside). 
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Fig.  4.10  Averaged  profiles  of  salinity  (southside). 
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Fig.  4.1 1  Averaged  profiles  of  salinity  (westside). 
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Fig.  4.12  Averaged  profiles  of  sigma-t  (northside). 
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Fig.  4.13  Averaged  profiles  of  sigma-t  (eastside). 
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Fig.  4.14  Averaged  profiles  of  sigma-t  (southside). 
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Fig.  4.15  Averaged  profiles  of  sigma-t  (westside). 
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Fig.  4.16  Vertical  section  of  temperature  for  circuit  1. 
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Fig.  4.17  Vertical  section  of  temperature  for  circuit  2. 
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Fig.  4.18  Vertical  section  of  temperature  for  circuit  3. 
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Fig.  4.19  Vertical  section  of  temperature  for  circuit  4. 
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Fig.  4.21  Vertical  section  of  temperature  for  circuit  6. 
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Fig.  4.22  Vertical  section  of  temperature  for  circuit  7. 
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Fig.  4.23  Horizontal  velocity  difference  plots  for  75-25m  levels  (circuit  2). 
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Fig.  4.24  Horizontal  velocity  difference  plots  for  75-25m  levels  (circuit  3). 
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Fig.  4.25  Horizontal  velocity  difference  plots  for  75-25m  levels  (circuit  5). 
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Fig.  4.26  Horizontal  velocity  difference  plots  for  75-25m  levels  (circuit  6). 


55 


Fig.  4.27  Horizontal  velocity  difference  plots  for  75-25m  levels  (circuit  7). 
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Fig.  4.28  Vertical  section  of  salinity  for  circuit  1. 
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Fig.  4.29  Vertical  section  of  salinity  for  circuit  2. 
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Fig.  4.30  Vertical  section  of  salinity  for  circuit  3. 
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Fig.  4.31  Vertical  section  of  salinity  for  circuit  4. 
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Fig.  4.32  Vertical  section  of  salinity  for  circuit  5. 
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Fig.  4.33  Vertical  section  of  salinity  for  circuit  6. 
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Fig.  4.34  Vertical  section  of  salinity  for  circuit  7. 
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Fig.  4.35  Staggered  profiles  of  temperature. 
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Fig.  4.36  Staggered  profiles  of  salinity. 
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Fig.  4.37  Staggered  profiles  of  sigma-t. 
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Fig.  4.38  Vertical  section  of  sigma-t  for  circuit  1. 
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Fig.  4.39  Vertical  section  of  sigma-t  for  circuit  2. 
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Fig.  4.40  Vertical  section  of  sigma-t  for  circuit  3. 
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Fig.  4.41  Vertical  section  of  sigma-t  for  circuit  4. 
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Fig.  4.42  Vertical  section  of  sigma-t  for  circuit  5. 
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Fig.  4.43  Vertical  section  of  sigma-t  for  circuit  6. 
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Fig.  4.44  Vertical  section  of  sigma-t  for  circuit  7. 
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Pig.  4.45  Average  profiles  of  the  static  stability  parameter  (E) 

for  the  northside. 
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Fig.  4.46  Average  profiles  of  the  static  stability  parameter  (E) 

for  the  eastside. 
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Fig.  4.47  Average  profiles  of  the  static  stability  parameter  (E) 

for  the  southside. 
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Fig.  4.48  Average  profiles  of  the  static  stability  parameter  (E) 

for  the  westside. 
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Fig.  4.49  Average  profiles  of  the  Richardson  number  (northside). 
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Fis.  4.50  Average  profiles  of  the  Richardson  number  (eastside). 
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Fig.  4.51  Average  profiles  of  the  Richardson  number  (southside). 
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Fig.  4.52  Average  profiles  of  the  Richardson  number  (westside). 
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Fig.  4.53  Staggered  profiles  of  the  static  stability  parameter  (eastside). 
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Fig.  4.54  Staggered  profiles  of  the  static  stability  parameter  (westside). 
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V.  SUMMARY  AND  CONCLUSIONS 


Water  column  stability  is  determined  by  the  interaction  of  dynamic  and  static 
forcing.  Examination  of  the  stability  and  the  driving  mechanisms  in  the  upper  100m  of 
the  measurement  area  was  accomplished  by  the  graphical  display  of  the  temperature, 
salinity,  sigma-t,  and  velocity  fields  as  well  as  the  static  stability  parameter  and  the 
Richardson  number.  The  analysis  of  these  displays  indicated  a  rapidly  changing 
environment  and  complex  horizontal  structure. 

Several  forcing  mechanisms  and  their  effect  on  stability  were  investigated  within 
the  context  of  the  analysis  fields  and  the  atmospheric  and  tidal  conditions  existing 
during  the  study  period.  These  mechanisms  were  upwelling,  currents,  tides,  bottom 
topography,  and  internal  waves.  All  of  the  mechanisms  were  found  to  be  a  factor  in 
the  stability  of  the  water  column  with  the  exception  of  upwelling.  The  California 
Current  was  found  to  dominate  with  an  apparent  absence  of  the  Davidson  Current. 

The  analyzed  fields  of  temperature,  salinity,  density,  and  velocity  exhibited 
considerable  spatial  and  temporal  variability.  This  variability  took  the  form  of  wave¬ 
like  patterns,  migrating  parcels,  fluctuating  isolines,  and  intensification/weakening  of 
gradients.  These  disturbances  could  be  attributed  to  the  forcing  mechanisms  of  bottom 
topography,  tidal  flow,  current  systems,  and  internal  gravity  waves.  The  shear,  as 
measured  between  25m  and  75m  was  calculated  to  be  at  least  5.0  cm  sec.  The  presence 
of  shear  provided  for  the  possibility  for  dynamic  instability.  Double  diffusive  processes 
were  indicated  to  be  a  factor  in  stability  by  the  small  reversals  in  the  density  field 
revealed  in  the  staggered  static  stability  profiles.  However,  despite  the  large  variability 
evident  in  the  temperature  and  salinity  fields,  the  sigma-t  field  did  not  exhibit 
comparable  variability.  This  indicates  that  most  of  the  temperature  and  salinity 
disturbances  were  density  compensating  and  thus  maintained  stability. 

The  stability  of  the  water  column  was  found  to  be  both  statically  and 
dynamically  stable  on  the  average  as  would  be  expected  considering  the  episodic  nature 
of  instability.  However,  small  patches  of  instability  were  identified  in  the  staggered 
profiles  of  static  stability  which  were  predominantly  located  in  the  25m  to  75m  region 
of  the  water  column.  This  region  of  the  water  column  is  also  where  the  greatest 
variability  in  the  thermohaline  and  velocity  fields  occured.  The  patches  of  instability 
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appeared  to  be  primarily  due  to  double  diffusive  processes  initiated  by  periodic 
intrusions  of  warmer,  less  saline  water. 
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